Stabilizer Codes for Continuous-variable Quantum Error Correction 
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A continuous- variable generalization of the discrete- variable stabilizer formalism for quantum error 
correcting codes is presented. This generalization is a step toward an independent understanding 
of continuous-variable quantum information. Our formalism yields all continuous-variable codes 
discovered to date and can be used to construct continuous-variable analogues of discrete- variable 
and classical codes. We use it to rederive the nine-mode code given by Lloyd and Slotine, and a five- 
mode code obtained by Braunstein. In addition, we construct a new continuous- variable code based 
on a code of Gottesman mapping three logical modes of information onto eight physical modes. 



One of the necessary techniques for the realization of 
quantum computing is a way to encode and manipulate 
quantum information in a fault-tolerant manner, so that 
small, unintended perturbations in the state of system 
have little or no effect on its large-scale dynamics. To en- 
able this, quantum error-correcting codes, and in particu- 
lar, stabilizer codes have been developed. For systems in 
n-dimensional state space, especially for qubits [n = 2), 
general techniques for constructing quantum codes and 
their associated logical operators are well known. 

Continuous variables are a promising new flavor of 
quantum information whose potential is still largely 
unexplored. To date, there have been only a few 
continuous-variable error-correcting codes constructed. 
Lloyd and Slotine have generalized Shor's 9-qubit code 
0, and Braunstein constructed a five-mode code, and 
conjectured about a general technique for converting 
discrete- variable quantum codes into continuous- variable 
quantum codes 0] . 

In this paper, we present such a general technique. 
Just as Gottesman Q and Calderbank et al. Q showed 
a correspondence between Pauli operators on n qubits 
and vectors in GF(2) 2n , we show an correspondence 
between operators in the Heisenberg-Weyl (generalized 
Pauli) group and vectors in M. 2n . By rewriting vec- 
tors in GF(2) 2n as vectors in R 2 ™, we have a general 
method of constructing continuous variable codes from 
discrete-variable codes. As an example, we can recon- 
struct the examples presented earlier, and present a new 
code based on a code of Gottesman which encodes three 
logical qubits in eight physical qubits and corrects one 
error. 

The Pauli group G = {I,X,Y,Z} and the n-qubit 
Pauli groups G n — G® n play a large role in quantum er- 
ror correction, because they are a basis for operators on n 
qubits, and any code that can correct errors in a subset of 
the Pauli group can then also correct an arbitrary linear 
combination of such errors The continuous- variable 
generalization of the Pauli group is the Heisenberg-Weyl 
group, or generalized Pauli group. Let q and p be two 
canonically conjugate observables, [q,p] = 1, with contin- 
uous spectra, for example the position and momentum of 
a one-dimensional system. The corresponding state space 



has as a basis the spectrum of q, namely the "position 
eigenvectors" \q), q £ E. The generalized Pauli group is 
the collection of all linear shifts in position and momen- 
tum. It is generated by shifts in p and q, X{t) = e mt v 
and Z(t) = e 1 ^ . Thus, an arbitrary Pauli operator has 
the form X(s)Z(t) for some real numbers s and t. The 
generalized Pauli group acting on n of these systems to- 
gether is then the collection of tensor products of these 
operators 



<X(s n )Z(t n ) : 8i,ti (1) 



A typical member of this group can expressed in terms of 
the canonical observables qi and pi as a unitary operator 
indexed by a vector v = (si, . . . , s n , tt, . . . , t„) G R 2 " 

U(v) = exp ^^(s^ + tiQiU , (2) 

with any two operators commuting up to a phase: 

U(v)U{v') = e i7ra, ( v ' v V(v')£/-(v), (3) 

where a>(v, v') = X^i^^i — s i^») ^ s the standard sym- 
plectic form on R 2n . This gives us an immediate corre- 
spondence between Q n and the symplectic vector space 
R 2 ™: The operator U(v) corresponds to the vector v, 
multiplication of operators corresponds to addition of 
vectors (up to phase), and scaling of vectors by real num- 
bers corresponds to taking real-number powers of opera- 
tors. 

Furthermore, we can use vectors in these symplectic 
spaces to build up abelian subgroups of the Pauli group: 



Suppose we have k vectors u*, . . . , u% G 



p 2/) 



such that 



<x>(iii, Uj) = for all i, j. Then we have a subgroup 



S = {C/(v) : v = ^OiUi, a. t £ 



(4) 



of Q n in one-to-one correspondence with the span of the 
vectors u^. This subgroup S is actually abelian, since any 
two operators U(v) and U(v') commute if w(v,v') = 0, 

and if v = X)i=i a i u i an d v ' = X«=i a i u n tnen 



w(v,v') = a t a' ? ujju^Uj 







(5) 
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The construction of stabilizer codes was given by 
Gottesman in Given a subgroup S of the n-qubit 
Pauli group such that (1) any two elements of S com- 
mute, and (2) —I is not in S, then the space of states 
stabilized by S consitutes a nontrivial code 

C(S) ={\ip): M \tp) = |V>> ,VM € S}. (6) 

The same fact is true in continuous dimensions: Any sub- 
group S of the continuous- variable Pauli group satisfying 
(1) and (2) stabilizes a nontrivial code. The case when 5* 
is a discrete subgroup was studied by Gottesman et al • 
They showed that such a discrete subgroup stabilizes a 
finite-dimensional subspace of the underlying state space, 
so that the codes created with such stabilizers cannot 
be used to encode continuous- variable quantum informa- 
tion. 

The stabilizers constructed in [j| correspond to sym- 
plectically integral lattices, namely the set of integer lin- 
ear combinations of 6 R 2n such that w(ui,u.,) is an 
integer. This criterion is too weak to allow for continu- 
ous subgroups, since such subgroups would allow for real 
linear combinations of the basis vectors. If, on the other 
hand, we take any set of k vectors u\, . . . , u k 6 R 2n sat- 
isfying the further restriction that uj(ui,Uj) = for all 
then as described in part II, the subgroup S corre- 
sponding to these vectors is an abelian subgroup of the 
continuous- variable Pauli group Q n , and thus a valid sta- 
bilizer. 

Since any operator which commutes with everything in 
S maps codewords to codewords, the group of such op- 
erators, namely the normalizer N(S) of S, has a natural 
interpretation as a collection of operators acting on the 
code space. In our symplectic language, an observable 
commuting with S corresponds to a vector v such that 
o;(v, Uj) = for all i, so that if the Ui span a subspace 
W C R 2 ™, then the space of operators commuting with S 
is the symplectic orthogonal W u = {v : cu(v, u^) = OVi}. 
It follows that dim = dimR 2 ™ - dim = 2n - k, so 
that the group N(S)/S of operators acting nontrivially 
on the code space corresponds to a subspace of dimension 
dimW u /W = dim -dimW = 2(n-k). Now we can 
construct a basis Xi, zi, . . . , x n _fc, z„_fc e R 2 ™ for /W 
such that w(xi,Xj) = uj{zi,%j) = and uj(x.i,Zj) = Sij. 
The group of operators corresponding to this basis satis- 
fies [U(sXi),U(tXj)} = [Z7(szO,^(*Zi)] = and 

U( SXl )U(t Zj ) = e^^Uitx^Uisxi), (7) 

and thus has the same group structure as the general- 
ized Pauli group G n -k, with the operators U(sXi) and 
£7(szj) acting as position and momentum shifts on the 
i-th logical mode, respectively. 

This means that if we can find k vectors u\, . . . ,Uk £ 
R 2 ™ such that u(ui,Ui) — 0, then we can choose any 
vector in the code space to represent the n — fc-mode state 
... 0), and use the n — k logical operators to transform 



this logical zero state into an encoding of any other state. 
In brief, k symplectically orthogonal vectors give rise to 
a continuous-variable code encoding n — k continuous- 
variable systems. 

Note that the states which are codewords stabilized 
by S have certain observable properties: Since they are 
stable (i.e. eigenvectors with eigenvalue 1) under the ac- 
tion of operators of the form U(xij) exp(^ i Sijpi + Ujqi), 
they must be eigenvectors of the observable ^ i SiPi + ttqi 
with eigenvalue zero. Thus, for each basis vector Uj = 
(s%j, . . . , s n j, tij, . . . , t n j) we have an observable mj = 
s ijPi +tij<h which is zero on the code space, and con- 
versely, the code space is the set of eigenvectors of these 
obsevables with eigenvalue zero. Thus, if we can measure 
these observables, which are just linear sums of position 
and momentum, then we can detect whether or not a 
given state is in the code stabilized by S. 

Based on the above, we have a procedure for construct- 
ing a continuous- variable quantum code: In order to cre- 
ate a code encoding n — k logical modes into n physical 
modes of information, we follow the following steps 

1. Choose k vectors Ui,...,Ufe G R 2 " such that 
w(uj,Uj) = for all 

2. Compute the standard logical zero state 

|0> = (E M)|0...0> (8) 

Mes 

= J dtf- dt k U(h ui) • • • U(t k u fe ) |0 . . . 0) (9) 

This state is obviously stable under the action of 
operators in S. 

3. Compute a set of logical Pauli operators. For these 
we can take any basis for the symplectic orthogonal 
W u of subspace W C R 2 ™ spanned by the vectors 
Ui, . . . ,u fc . 

Note that all three of these steps are nearly trivial if we 
already have a discrete- variable stabilizer code encoding 
n — k logical qubits into n physical qubits. In the first 
step, consider the binary vectors defining the discrete- 
variable code as real vectors with entries in {0, 1} C R, 
ensuring that w(uj,Uj) = 0, (often possible just by 
changing some entries from 1 to —1). The second step is 
a trivial computation in either case. In the third step, the 
binary vectors defining logical operators on the discrete- 
variable code can again be considered as real vectors, 
and appropriately modified so that they commute with 
the continuous- variable stabilizer. 

The smallest possible examples of error-correcting 
codes are the three-mode position and momentum codes 
mapping \q) to 

\q) = \qqq) , \q) = J d*l dt 2 dt 3 e **fr+*>+t>)9 \ tl ,t 2 , t 3 ) , 

(10) 
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and correcting a single shift in q or p, respectively. Ob- 
serve that the position code is invariant under the oper- 
ators Z(t) Z(-t) I and I Z(t) Z(-t): 

(Z(t) Z(-t) I) \qqq) = e "'V^~^ \qqq) = \qqq) . (11) 

Conversely, it is the code stabilized by the group S = 
(Z(t)Z(-t)I,I Z{t)Z(-t)) generated by these opera- 
tors, which corresponds to the generator matrix 



A = 








-1 

1 -1 



(12) 



This is a valid code, since we only have two genera- 
tors for the stabilizer, corresponding to the rows Ui = 
(0,0,0,1,-1,0) and u 2 = (0,0,0,0,1,-1) of A, and 
w(ui,U2) = 0. The logical zero generated in the stan- 
dard way is also what we expect it to be, namely 



|0) = / dh dt 2 e 



iir(ti—t-2)q 



000) 



d(t 1 -t 2 )S(t 1 -t 2 ) |000) = 000) 



(13) 
(14) 



The standard logical position shift, X(t) = X(t)(&X(t)® 
X(t), is still valid here, since its vector x = (1, 1, 1, 0, 0, 0) 
satisfies w(x, ui) = w(x, u 2 ) =0 for the basis vectors Ui, 
U2 of the stabilizer. This means that an arbitrary state 
\q) is encoded as X(q) |0) = (X (t) ® X (t) ® X (t)) |000) = 
\qqq), which is the standard definition of the three- mode 
position code. A similar construction, reversing the roles 
of X and Z , gives us a three- mode stabilizer code protect- 
ing against momentum shifts, and encoding an arbitrary 
position eigenvector as shown above. 

These two codes can be concatenated in the usual way, 
by first encoding one qubit into three with the position 
code, then encoding each of the three bits with the mo- 
mentum code. This gives Lloyd and Slotine's generaliza- 
tion of Shor's code, which encodes the state |g) as 

\q)= [ e"«( tl+t2+t3 > Mi,ti,t 2 ,t2,t2,t3,t3,t3> ■ 



(15) 

And the logical operators acting on the code, as in the 
case of the Shor code, are X(t) = {X{t)) m and Z(t) = 
(Z(t)) m . 

Braunstein has presented a five-mode code by showing 
a network which performs encoding this can also be 
described in this stabilizer formalism. By examining the 
logical zero state given by Braunstein, we find that its 
stabilizer corresponds to the generator matrix 



f 1 








-1 


-1 








1 


-1 


°\ 





1 








1 


1 








1 











1 


1 





1 











1 


V° 

















-1 


1 


-1 





A = 



so that we have the encoded zero state 

|0) = / dt 1 dt 2 dt 3 e i ^ lt2 \t 3 ,t 2 MM -*3,*2 



(17) 



*3>. 



For the logical Pauli operators, we take X(t) = Z(t) <g> 
X(t) ® X(t) ® I ® / (in order to replicate Braunstein's 
code), and take Z(t) = (Ft)® 5 X(i)(F)® 5 = X(t)®Z(t)® 
Z(t) ® I <g> I to be the conjugate of X(t) by a fourier 
transform acting on each mode. These logical operators 
correspond to the vectors x = (0, 1, 1, 0, 0, 1, 0, 0, 0, 0) and 
z = (1, 0, 0, 0, 0, 0, 1, 1, 0, 0), and thus commute with the 
stabilizer and form logical Pauli operators since w(x, z) = 

2 — 1 = 1. Acting with X(x) on the logical zero state, 
we find that the encoded version of the state \x) is the 
logical state 

\0 = (Z(q)®X(q)®X(q)®l2>r)\Q) (19) 
dU e **(*i*a+*»ff) \t 3 , ij, + q,tx + q, ti-t 3 ,t 2 - t 3 ) , 

(20) 

which is the same encoded state produced by the network 
described by Braunstein in 0. Thus, this stabilizer is 
another way of describing Braunstein's code, with the 
added advantage of providing the logical operations X{t) 
and Z{t). 

Finally, to show that this formalism can be used to 
produce new continuous-variable codes, we can general- 
ize a more interesting code of Gottesman, which encodes 

3 logical qubits in 8 physical qubits. The generator ma- 
trix corresponding to this stabilizer, shown in Table [I] is 
identical (except for signs) to the real generator matrix 
for the continuous-variable code. The logical Pauli op- 
erators follow similarly; by adding sign changes where 
necessary to the logical operations on the discrete code, 
we obtain the six operators acting as logical Pauli oper- 
ators on the encoded codewords, listed in Tabled Thus 
we can write the encoding of any three-mode position 
eigenstate in terms of these logical operators acting on 
the logical zero state: 



A = 



/llllllll 
00000000 
10 110 10 
10 10 10 1 

V 1 1 1 1 



000 \ 

111 1 -1 -1 -1 -1 ' 

0-1-1 1 1 

1-1 1-1 

10-1 1 -1 / 



(18) 



Xi(t) = X(t) <g> X(-t) ® i" «> / (8 1 ® Z(t) (8! I ® Z(-t) 

X 2 (t) = X(t)®I®X(-t)®Z(-t)®I®I®Z(t)®I 

X 3 (t) = X(t)®I®I®Z<g,X(t)®Z(t)®I®I 

Zi(t) = I®Z(-t)®I<8)Z(t)®I<8)Z(-t)®I<8)Z(t) 

Z 2 (t) = I®I®Z(-t)®Z{-t)®I®I®Z(t)®Z{t) 

Z 3 (t) = I ® I ® i" ® I ® Z{t) ® Z(t) ® Z(-t) ® Z(-tjjl6) 



TABLE I: Generator matrix and logical operators for the 
eight-mode code. 
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1 000) = / dU \tx,tx + t 3 +t A + t 5 ,t 1 +t 5 ,t 1 +t 3 + t 4 , h+t 3 + t 6 ,t 1 + t 4 , h + t 3 , h+t 4 + t 5 ) (21) 
mm) = (X 1 (q 1 )X 2 (q 2 )X 3 (q 3 ))\m) (22) 
= J dti e ^(9i*5+ 9 2t4+g3t3)| gi +q2 + q3 + tljtl + t3 + t 4 + h - qi,h +t 5 - q 2 , 

h+t 3 +t 4 + q 3 , ti+h + t b M+ U, h+t 3 ,t 1 +t 4 + t 5 ). (23) 

I 



so that we have a representation of the entire three-mode 
state space as superpositions of these eight-mode states. 
This is the first example of a continuous-variable error- 
correcting code which encodes more than 1 logical mode 
of information. 

Our generalization of Gottesman's stabilizer formalism 
provides an explicit algorithm for creating continuous- 
variable quantum codes out of discrete-variable quan- 
tum codes. In particular, it brings the power of the 
classical theory of error correcting codes to bear on 
continuous-variable problems, since any code created 
with the Calderbank-Shor-Steane construction, for ex- 
ample, can be translated into continuous-variable codes 
as in the above example. 

The creation of a large class of error-correcting codes 
for continuous- variable systems indicates several possibil- 
ities. First is that a high-level symmetry exists between 
the discrete- and continuous- variable theories of quantum 
error correction and fault-tolerant computation. Already 
in this vein, Bartlett et al. recently showed in that 
a version of the Gottesman-Knill theorem holds for con- 
tinuous variable systems, which describes how a large 
class of Hamiltonians in continuous dimension can be ef- 
fectively simulated on a classical computer, just as in 
discrete systems. Second, the combination of the codes 
as in 0, encoding discrete systems into continuous sys- 
tems, with the true continuous codes described above, 
could make possible much more interesting codes which 
combine systems of different dimensions (Braunstein con- 
jectured about these in [4j). These could, for example, 
be constructed as the codes stabilized by a subgroup of 
the continuous Pauli group which is discrete in some di- 
mensions and continuous in others. 

In practice, the approximation of position eigenstates 



by finitely squeezed states will hinder error-correction 
using continuous-variable codes. Braunstein noted in 
4] that this squeezing must make vacuum noise small 
with respect to both the expected size of errors and rel- 
evant length scales the states to be encoded. Gottesman 
et al. computed the effect of such dissipative error on 
their discrete- variable encodings in Q, but with tech- 
niques that do not immediately generalize to continu- 
ous codes. Similar bounds on the fidelity of continuous- 
variable codes with finite squeezing would allow quan- 
titative measures to be set on the absolute capacity of 
a continuous-variable channel, and could be the key to a 
theory of fault-tolerant continuous- variable computation. 

We would like to thank Olivier Pfister for his advice 
and guidance, and Harold Ward and James Rovnyak for 
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